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1 Introduction 



The decay channel Z — > WX was suggested in the early days of the Standard Model 
(SM) as a possible source of the PU-bosons [|IJ. A number of detailed studies arrived at 
the conclusion that the branching ratio of this decay mode is extremely small. There 
are several factors which contribute to this strong suppression. First, it is a higher order 
electroweak process; second, it is suppressed due to the small difference in masses of the Z 
and W bosons; and third, it suffers from a destructive interference among contributions of 
different diagrams. On the other hand it was also pointed out that such a decay process 
would lead to an exceptional signature (in the case when W boson decays leptonically). 
Namely, there will be one charged lepton with the transverse energy E± ~ 40 GeV 
accompanied by the same amount of missing energy. If the W decays hadronically it will 
be more difficult to isolate the signal from the background. 




Figure 1: Contributions to the decay Z — > Wud in the Born approximation 

An interesting feature of this process is that it involves the three gauge boson coupling 
ZWW (see fig. la). This coupling gives the largest contribution to the decay rate. Un- 
fortunately, there are also large gauge cancellations from the interference of the diagram 
with three gauge boson coupling with the other diagrams involved. In the SM the cou- 
pling ZWW is fixed. It is likely that effects of physics beyond the SM can modify this 
coupling and lead to the deviations from its SM value. In the context of the Z — > WX 
decay channel this possibility has been analyzed in Ref. @. It has been concluded that 
in a special case the introduction of anomalous couplings the rate can be increased by up 
to one order of magnitude. 

Let us present some numbers from Ref. concerning the number of events one can expect 
per 10 7 Z events in the Standard Model: 

N(Z -> Wlv -> lulu) = 0.19, 

N(Z -> Wlv -> qqlu) = 0.38, 
N(Z -> Wqq -> luqq) = 0.35, 
N(Z -> Wqq -> qqqq) = 0.70. 
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Hence to observe the decay Z — > WX within the SM requires more than 10 8 Z events. 
The largest sample of Z events has been collected by the four experiments at LEP1. 
Though the number of Z's is huge (~ 16 • 10 6 ) it is still not sufficient to observe the decay 
Z — > WX. On the other hand, the designed luminosity of a future e + e~ collider will 
allow to produce 10 7 Z bosons per day, if the machine operates on the Z resonance. Such 
experiments are indeed being planned at the JLC and will allow a study of the rare decay 
channels of Z and in particular of Z — > WX. 

In Ref. [|J the QCD corrections to this decay mode were discussed. The authors of 
that paper mainly concentrated on the non-perturbative and logarithmic perturbative 
corrections. They concluded that the considered corrections are small and can not 
strongly shift the lowest order branching ratio. Nevertheless, pure perturbative corrections 
0(a s (Mz — Mw)) remained unknown. 

In this paper we present a simple estimate of these corrections, based on the expansion 
of the rate with respect to the small ratio (Mz — Mw)/Mz- In the next section we 
demonstrate the basic idea of our approach with the example of the lowest order decay 
rate. Then we discuss the QCD corrections; our results are summarized in the Conclusions. 

2 Lowest order decay rate 

Let us describe the idea of the calculation with the example of the Born approximation. 
We consider Mz and Mw as large parameters (compared to Mz — Mw)- Then the quarks 
in the final state will have relatively small energies and momenta. Hence we can try to 
expand the amplitude for Z — > Wud in the quantities of the order of (Mz — My/)/Mz- 

We adopt the following notation: p and k are momenta of Z and W, p\ and P2 are 
momenta of the quarks in the final state. 

We begin with the diagram of Fig. lb where the H^-boson is emitted from a quark line. 
The virtual quark is far off-shell. Contracting its propagator to a point and neglecting 
momenta of the final quarks we obtain the following expression for the amplitude: 



Here L — (1 — 7s)/2 and g q _ describes Z coupling to the left-handed quarks, g q _ = 1%/sw 
Sw denotes the sine of the weak mixing angle. We note that according to our approxima- 
tion we put its cosine equal to 1. 

Using Dirac algebra identities we can rewrite this as 



M 2 



V2s w M^ 



u(p 1 )Ov(p 2 ) 



O = e z pe w L. 



(1) 



M 2 




T (C li + A lt )u( Pl )'fLv(p 2 ) 



(2) 



VlswMi 



w 




(3) 



e w ■ P = £z ■ P = 0, 



(4) 
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which is exactly valid if Myy = Mz- 

The interference of the two pieces of the amplitudes, proportional to and A^, vanishes 
after the sum over W and Z polarizations. Note also that the amplitude is described 
by the same equation, only the sign of C M is opposite. 

The amplitude Aii involves a three gauge boson coupling. Taking there the limit P1+P2 — * 
we arrive at the approximate formula: 

ie 2 y/2 

Ml = 2 ti#2 ( e w ■ ez)Pnu(Pihv Lv (P2)- ( 5 ) 

To calculate the leading asymptotics of the Born width we integrate over the light quark 
phase space and then over W three-momenta. Because of the approximate equality 
Mz ~ M\y the W boson is always non-relativistic. We approximate its phase space 
element by 

* = % ( 6) 

2E k 2M W 2 ' M| v ; 

By counting the powers of (Mz — My/) we see that the leading order behavior of the 
tree-level width is given by (1 — x) 5 . 

The decay width for Z — > Wud can be written as 

1 2 1 d ^1 I „u „d \ 



Y(Z -> Wud) = T 



n u2 _|_ n d 1 ( n u _ n a 

9 - Z 9 - H x (x) + — H 2 (x) + g u _g d _H 3 (x) - ™= 9 -=± H 4 (x) 

l S w S]y 



In the above equation we denote 



(7) 

N c M z a 2 



167TS^ 



and Hi(x) are the contributions of various diagrams to the width. Namely, Hi(x) is the 
contribution of the squares of the graphs lb or lc, H2(x) is the contribution of the graph 
la squared, H%(x) is lb — lc interference and H±(x) is determined by la — 16 and la — lc 
interferences. 

The leading asymptotics of these functions are 



Hi(x) = (1 - x) 5 hr(l + C F ^SA; (9) 

1 20 2 30 3 60 4 30 
The results of a complete Born level calculation can be found in 0. Let us note that since 
we use the values = — g d _ = ^7 we obtain at Born level complete cancellation between 
H 2 and H^, potentially the largest contributions to the width. This is the destructive 
interference which makes the rate additionally suppressed. 

Our aim in this paper is to discuss the QCD corrections to this decay width which are 
parametrized by the coefficients Si. 
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3 QCD corrections 



First, consider QCD corrections to the square of the amplitude M.\. Since the gluons 
couple to the external quarks the only QCD correction to |.Mi| 2 will be the correction to 
the correlator of two V — A currents, which can be taken from VF-decay. This has been 
calculated for the first time in Ref. ||. The ratio of the one-loop correction to the Born 
contribution is[] a s /n. 

Next, we analyze the QCD corrections to the remaining amplitudes. As an example let 
us take the graph were P^-boson is emitted from the down quark line (M.2). We discuss 
first the radiation of real gluons. 

There are three graphs describing real gluon emission. In two of them the gluon is 
emitted from the external fermion line and in the third one from the fermion line with 
high (C(M^)) virtuality. The energy of the gluon is restricted due to the phase space 
suppression, therefore the third graph will have an additional suppression factor 0((Mz — 
Mw) 2 /Mz) relative to the Born one. It is therefore of no interest for us. 

The leading contribution of the two remaining graphs can be found by contracting the 
virtual fermion propagator carrying the momentum 0(M\y,z) to a point. The effective 
vertex which appears as the result of such contraction is equivalent to the effective vertex 
presented in the eq. (|]). Unfortunately, the transformation which has been used in the 
transition from eq. (JI|) to eq. (|2]) can not be applied for the graphs which describe radiative 
corrections. The reason for this is the use of dimensional regularization and the presence 
of infra-red divergences in the graphs with emission of real gluons. This problem makes 
it necessary to perform an "honest" calculation of the 0(a s ) correction to the production 
of hadrons by the operator O defined in eq. ([I]). 

Next, it is necessary to analyze virtual corrections. There are four of them. Note that 
there are no self energy corrections to the external quark lines because they are represented 
by no-scale diagrams which vanish in dimensional regularization. If we perform on-shell 
renormalization for the external quark lines this also means that there is no wave-function 
renormalization at all. Therefore the sum of virtual and real corrections must be finite 
on its own. 

There is a simple way to rearrange virtual contributions to make them more transparent. 
First let us note that in the vertex corrections to Z —>■ u*u or u* — > Wd we can simply 
put the momenta of the "soft" quark equal to zero; this will not introduce any new 
divergences. The same also applies to the self energy correction to the virtual fermion 
line. 

The remaining contribution is due to the box diagram. Here the situation is slightly more 
delicate — we can not put the momenta of the light quarks equal to zero because this leads 
to additional divergences. The solution is to apply the ideas of asymptotic expansions 
@. In that approach the contribution of the box graph can be obtained as a sum of two 
pieces: in the first one the large momentum flows through the fermion propagator which 
connects the Z and W vertices. To get the leading contribution we must contract it to 
the point. The second piece is the expansion of the whole box in Taylor series in the small 
external momenta. In both of these pieces one gets additional divergences (ultraviolet in 

3 This result is true not only for the leading asymptotics in Mz — M\y- 
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the "contracted" part and infrared in the part originating from the Taylor series). These 
divergences are canceled in the sum. However, it is very useful to separate these pieces 
because of the following observation: The "contracted" part of the box taken together 
with the real emission graphs will give a complete 0(a s ) correction to the production 
of hadrons by an operator O in eq. (jl]). Therefore, this sum should be ultraviolet and 
infrared finite. For obvious reasons we will call this contribution "soft". 

The remaining ("hard") virtual corrections, i.e. corrections to the Z —>■ u*u vertex, u* — > 
Wd vertex, self energy correction to the u* propagator, and the box graph, all taken at 
zero momenta of external fermions, must also be finite. Note, however, that individual 
pieces are divergent and hence one has to perform the Dirac algebra in D-dimensions. 
The final expression for the hard correction to the amplitude which can be obtained in 
this way is remarkably simple: 



To get an expression for Ai^ ard it is sufficient to reverse the sign in front of C M in eq. fllPl). 
One gets 

M >" d = Tfik^ 2C " + 7 ^«7^M. Hi) 

We remind that there is no interference between C M and structures, hence if we know 
their relative contributions to the Born width it is a trivial exercise to find a value of the 
QCD correction. 

Let us demonstrate how this works by considering the QCD corrections to the interference 
of the diagrams Aii and Consider first the correction due to the soft part of the 

diagram The operator which "produces" final hadronic state can be written as 

e w pe z L = -p [e w ■ e z + -[ew, ^z]j L. (12) 

Here we have used Eq. ([|). Note that because of the sum over polarizations of the Z and 
W bosons there is no interference between the first and the second structure in the above 
equation. The amplitude .Mi is always proportional to the product ew ■ e Z - Therefore in 
the soft correction only the first structure in eq. fll2|) contributes. The correction to it is 
again just the correction to the W decay width. As a result we find 

S? ft = | (13) 

Now the "hard" part of the correction comes from the interference of the C M structure of 
eq. (|TT|) and eq. (||]). Comparing coefficients and signs of the "hard" correction with the 
corresponding Born one we conclude that the hard correction to H^(x) is 

S hard = _ 1 (M) 



and 



5 4 = 5l oft + 5 h 4 ard = 1 (15) 
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Clearly, the contribution of the diagram with the down-type quarks coupling to W^-boson 
is the same up to trivial redefinitions of the coupling constant. 

In a similar manner we obtain the QCD correction to the square of the amplitudes }A 2 
or jVf 3 and for their interference. The hard correction can be immediately determined 
from M.2% d - F° r the soft part, however, one needs a full calculation. This is by no means 
difficult and technically is equivalent to the calculation of the QCD correction to the e + e~ 
annihilation to massless quarks. All together this gives 

Si = ~ 4 (16) 



We find also 



(17) 



Let us summarize the results of the QCD corrections to the quantities Hi(x): 



^hard . 
gsoft 



Shard 
°2 



_8 

3' 

25 rso/t 
12' °2 



0. 

3 
•1' 



Shard 
°3 

ssoft 

°3 



-6, 



Shard 
°4 



19 ssoft 
T> °4 



_i 

2' 

3 
1' 



The division of the QCD corrections into the soft and hard parts helps us to determine 
the proper energy scale of the strong coupling constant. From the explicit calculations 
which we have described it is clear that the characteristic virtualities of the gluons are 
[Mz — Mw)/2 and Mz in the soft and hard parts, respectively. Therefore, instead of 
eq. (0) we use 



HAx) 



1-xYh 



5a(o) 



1 + C 



V 



a s (MJ 



r 



71 



zhard 



+ c t 



M z -M w 
2 



rsoft 



71 



(19) 
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Among the four correction factors 5i = S^ ard + 51° two can be compared with previously 
known results. 62 is the well known QCD correction to (axial) vector currents. 5i, on the 
other hand, is the first term of the expansion of the exact result for the square of the 
diagram lb (or lc) which was recently found in a rather unusual way 0. The complete 
mixed electroweak and QCD corrections to the Z boson decay width were computed in 
the limit of very small and very large W boson mass. Their difference gives precisely 
the QCD corrections to the emission of real W. From the several known terms of both 
expansions the general (very simple) formula for the coefficients was guessed; the exact 
sum of this series is, however, rather complicated: 



il6+lc 



a s _ gl 2 + g d2 



Cp — r 

7T 



7 1 3 
1 — In x H 

4 2 32a; 



9x. 2 /7 1 

— mx + x — I — In a; 
8 U 2 



.^l + ^ S (x)--S 
32 3 v ; 3 



(20) 



Here 



S{x) 



4 2 
— (x 2 + Ax + l) Li 4 (x) + — (7x 2 + 16x + 7) Li 2 i 



x 



7 



-f} (!-*•) hp 

+ ln(x 



137 



40 

x 



4 6 
- (x 2 + 4x + l) Li 3 (x) + 



X/ 



1-x 2 ) Li 2 (x) - 18- 



+ ln 2 (x) 



10 

a; 
23 



1 3 23 4 

~2 (x 2 + 4x + l) Li 2 (x) + — (l - x 2 ) ln(l - x) + — + 



x 



(21) 



where 1^2,3,4 (x) are polylogarithms ||. An expansion of this formula around x = 1 
confirms our result for 5j. One can expect that the exact result for the correction to 
the interference of diagrams lb and lc would be even more complicated and its full 
calculation is a rather daunting task. Therefore it is appropriate to use an expansion 
method described in this section. 



4 Conclusion 



We have presented a simple estimate of the 0(a s ) correction to the decay rate Z — > 
W~ud, W + ud. Our approach is based on the expansion of the event rate in powers 
of (Mz — M\y)/Mz- In principle, this expansion parameter is not very small and the 
leading term in the expansion is not expected to approximate the exact result with high 
accuracy. However, it gives us reasonable estimates of the lowest order and one-loop QCD 
corrections to the decay width almost without effort. Another important point is that 
the accuracy of this approximation can be improved by calculating further terms in the 
expansion, should need arise. Such a calculation is by no means impossible. 



Combining all QCD corrections presented in eq. (18) we arrive at the following correction 
to the decay rate: 



T(Z -> Wud) = T Born (Z -> Wud) 



4a^Mf) 

3 7T 



a. 



+ 



M z -M w 
2 



:22i 



7T 



where we have used Cf = 4/3. 

For the values of the strong coupling constant we use a s (M|) = 0.12 and a s 

= (0.25.. 0.30). Numerically this gives the correction to the decay rate Z - 
order of 13-15 percent. 



M z -M w 
2 



Wud of the 



Let us have a closer look at the accuracy of our formulas. We can compare our approx- 
imate results (0^|) for the lowest order event rate with the results based on a complete 
calculation [[J. We conclude that our results based on the leading asymptotics to the 
event rate give a 30 percent accuracy. Note also that the accuracy with which the ra- 
tios Hi/Hj are reproduced is even better and amounts to approximately 10%. We can 
therefore expect that the QCD corrections derived in this paper represent the complete 
unknown QCD corrections to the decay rate Z — > Wud with the accuracy 10-20 percent. 

Our final result is that the QCD corrections increase the decay rate of the process Z — > 
Wud by approximately 14 ± 2 percent. 
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